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ABSTRACT 


We present the Deep Synoptic Array (DSA-110) discovery and interferometric localization of the so 
far non-repeating FRB 20220319D. The FRB originates in a young, rapidly star-forming barred spiral 
galaxy, IRAS 020444-7048, at just 50 Mpc. Although the NE2001 and YMW16 models for the Galactic 
interstellar-medium (ISM) contribution to the DM of FRB 20220319D exceed its total observed DM, 
we show that uncertainties in these models accommodate an extragalactic origin for the burst. We 
derive a conservative upper limit on the DM contributed by the circumgalactic medium (CGM) of 
the Milky Way: the limit is either 28.7 pc cm"? and 47.3 pc cm" ?, depending on which of two pulsars 
nearby on the sky to FRB 20220319D is used to estimate the ISM DM. These limits both imply that 
the total Galactic CGM mass is < 10!! Mọ, and that the baryonic mass of the Milky Way is € 60% of 
the cosmological average given the total halo mass. More stringent albeit less conservative constraints 
are possible when the DMs of pulsars in the distant globular cluster M53 are additionally considered. 
Although our constraints are sensitive to possible anisotropy in the CGM and to the assumed form of 
the radial-density profile, they are not subject to uncertainties in the chemical and thermal properties of 
the CGM. Our results strongly support scenarios commonly predicted by galaxy-formation simulations 
wherein feedback processes expel baryonic matter from the halos of galaxies like the Milky Way. 
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1. INTRODUCTION 


Galaxies like the Milky Way are embedded in a multi- 
phase (~ 104 — 10" K), highly ionized (hydrogen neu- 
tral fractions « 0.0196), kinematically complex, spa- 
tially clumpy and anisotropic circumgalactic medium 
(CGM; Tumlinson et al. 2017). Likely extending beyond 
dark-matter halo virial radii ryi,, the CGM may repre- 
sent the dominant baryon component by mass within 
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halos. The physical properties of the CGM of exter- 
nal galaxies, including density, temperature, kinematic 
structure, and chemical composition have long been 
probed by absorption-line measurements towards back- 
ground objects. More recently, detections of thermal 
bremsstrahlung (e.g., Li et al. 2018) and the Sunyaev- 
Zeldovich effect around nearby galaxies (e.g., Bregman 
et al. 2022) provide independent constraints on the 
CGM density and temperature structure, potentially 
confirming the presence of extended coronae around 
galaxies. Evidence for extended gas reservoirs associ- 
ated with nearby galaxies, either in the CGM or the 
intra-group medium, has also been observed in the dis- 
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persion measures (DMs) of background fast radio bursts 
(FRBs; Connor & Ravi 2022; Wu & McQuinn 2022). 
Galaxy-formation simulations highlight the dependence 
of CGM properties on feedback from stellar winds, su- 
pernovae, and AGN, as well as on the galaxy merger his- 
tories (e.g., Wijers et al. 2020; Hafen et al. 2020; Zheng 
et al. 2020; Fielding et al. 2020; Appleby et al. 2021; 
Ramesh et al. 2022). The simulations predict varying 
degrees of anisotropy in the CGM of individual galax- 
ies, and different total CGM masses for galaxies with 
comparable global characteristics. A ubiquitous predic- 
tion of simulations is that feedback processes determine 
the integrated CGM mass, Mcam. In general, CGM 
baryon fractions foam = Moonta < 0.5 are predicted, 
where Mo is the total (dark matter and baryonic) halo 
mass. 

A halo of ~ 109 K gas at densities < 107? cm"? sur- 
rounding the Milky Way was proposed by Spitzer (1956) 
as a solution to the problem of confining distant cold- 
gas clouds at high Galactic latitudes. Tentative evi- 
dence for such coronal gas had previously been found 
in observations of the diffuse radio-synchrotron back- 
ground (e.g., Baldwin 1956). A theoretical basis for 
the origins of such a hot CGM around galaxies was de- 
veloped in the 1970s (Rees & Ostriker 1977; White & 
Rees 1978), with the gas infalling onto dark-matter ha- 
los heated to virial temperatures of > 109 K, possibly 
in shocks at the halo virial radii (see also Birnboim & 
Dekel 2003). It would, however, be nearly two decades 
before further observational evidence for a hot Milky 
Way CGM was identified using observations of shadows 
in the soft X-ray background from high-latitude (Snow- 
den et al. 1991) and high-velocity (Herbstmeier et al. 
1995) cold clouds, and in the spectral decomposition of 
the background into absorbed and unabsorbed compo- 
nents (Garmire et al. 1992). Today the hot CGM of 
the Milky Way is best traced through X-ray observa- 
tions of OVII and OVIII diffuse emission, and absorp- 
tion towards bright background sources. Although high 
Galactic latitude data are well fit by a spheroidal com- 
ponent with a mass of a few x10!°Mo (e.g., Bregman 
& Lloyd-Davies 2007; Gupta et al. 2012), more recent 
models include a substantially less massive (~ 0.1%), 
but denser disk-like component to explain the line ob- 
servations (Li & Bregman 2017; Nakashima et al. 2018; 
Kaaret et al. 2020; Ueda et al. 2022). Structure in the 
hot Milky Way CGM is also observed in the form of 
the Fermi (Su et al. 2010), Parkes (Carretti et al. 2013) 
and eROSITA (Predehl et al. 2020) bubbles. The neu- 
tral and warm (X 10? K) components of the Milky Way 
CGM, as primarily traced through observations of high- 
velocity clouds, are likely sub-dominant in mass to the 


hot component (few x105 Mo; Putman et al. 2012), but 
highlight the complexities of the disk-halo interface (e.g., 
Koo et al. 1992; Werk et al. 2019; Clark et al. 2022). 

'The mass of the Galactic CGM is uncertain. Here we 
assume Mio, = 1.5 x 102 Me (following Prochaska & 
Zheng 2019) and a total Galactic mass in stars and cold 
gas of 6 x 10^ M (Draine 2011). Setting foam = 1 
and (2/8 = 0.188 (Planck Collaboration et al. 2020) 
implies Meam = 2.82 x 101! Mo. Recent models for the 
distribution of OVII and OVIII emission and absorp- 
tion result in Mcaw values between 3 — 4 x 10!9 Mo 
(Li & Bregman 2017), 5.5 — 8.6 x 1019 Mc (Kaaret et 
al. 2020), and ~ 1.2 x 10 M (Yamasaki & Totani 
2020). Faerman et al. (2017) argue that a portion of 
the hot CGM is also traced by OVI absorption, as is 
observed in external galaxies (Tumlinson et al. 2011), 
and find Moam = 1.2 x 10! Mo. Systematic effects 
that make the estimation of Mcaw difficult include the 
model uncertainties in foregrounds and in the spatial 
distribution and clumpiness of the gas, and the poorly 
constrained gas metallicity. An alternate measurement 
of Mcam was found by Salem et al. (2015) by model- 
ing the ram-pressure stripping of the Large Magellanic 
Cloud (LMC); a low value of (2.7 + 1.4) x 101 Mc was 
found. 

'The DMs of FRBs and pulsars within the Milky Way 
halo offer a unique probe of the content of the Galactic 
CGM (e.g., Anderson & Bregman 2010; Platts et al. 
2020). FRBs are transient radio emissions from distant 
extragalactic sources that, like the emissions from radio 
pulsars, are so short in duration that the dispersion in 
intervening plasma is evident as an arrival-time delay at 
lower frequencies: 
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where qe is the electron charge, Mme is the electron mass, 
v is the frequency, D is the source distance, and ne(l) is 
the electron number density along the sightline. Under 
physical conditions of CGM gas, DMs quantify the col- 
umn densities along the source sightlines. In all cases, 
however, the observed DMs of FRBs, DMops, are com- 
posed of a selection of additive components: 


DMobs = DMism + DMcam + DMigm + DMhost, (2) 


where in this formulation DMyjsy includes the warm ion- 
ized medium (WIM) in the Milky Way disk (e.g., Cordes 
& Lazio 2002; Yao et al. 2017), DMcay is the DM con- 
tributed by the Galactic CGM, DMmpos, is the DM as- 
sociated with the FRB host-galaxy interstellar medium 
(ISM) and their halos, and DMyem includes contribu- 
tions from the intergalactic medium (IGM) and any in- 
tervening galaxy systems. For pulsars within the Milky 
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Way halo, such as in the Magellanic Clouds and in sev- 
eral globular clusters, DMjqq = 0pccm-^?. These ex- 
pressions apply to the redshift ~ 0 regime relevant to 
this paper. If the components of DM op; can be accu- 
rately modeled or bounded, constraints can be placed 
on DMcem. Modeling of the DMs of pulsars in the 
LMC led Anderson & Bregman (2010) to find a very low 
Moaw ~ 1.5x 101? M assuming a Navarro et al. (1997) 
profile (an NFW profile); it is likely that an overly large 
DMrzsyw was assumed, and the NFW profile is no longer 
favored for the Galactic CGM. Although most FRBs are 
too distant to allow for accurate modeling of the compo- 
nents of DMops, FRBs from nearby galaxies (Bhardwaj 
et al. 2021; Kirsten et al. 2022) have begun to be used to 
place constraints on DMcem, assuming specific models 
for DMjsy, that are in tension with some models for the 
Milky Way CGM. 

In this paper we present observations with the 
Deep Synoptic Array (DSA-110) of a new FRB 
(FRB 20220319D) that motivates stringent new con- 
straints on the mass of the Milky Way CGM. Remark- 
ably, FRB 20220319D was observed to have DMops < 
DMrsy for leading models of the warm-gas distribution 
in the Galactic disk, despite the unambiguous associ- 
ation we find with a galaxy at a distance of 50 Mpc. 
We detail observations of the FRB in 82. We con- 
sider the association with its host galaxy in 83, wherein 
we show that the DM of FRB 20220319D is consis- 
tent with an extragalactic origin given uncertainties in 
models for DMygm, using observations of the DMs of 
pulsars in globular clusters with accurate parallax dis- 
tances. We present an analysis of the host environ- 
ment of FRB 20220319D in 84. We then place conser- 
vative constraints on DMceam, and hence Moem, using 
FRB 20220319D and the distant high-latitude globular 
cluster M53, in 85. We discuss the implications of our 
results in 86, and conclude in 87. 


2. DSA-110 OBSERVATION OF FRB 20220319D 


The DSA-110! is a radio interferometer hosted at 
the Owens Valley Radio Observatory (OVRO), pur- 
pose built for the discovery and arcsecond-localization 
of FRBs. A full description of the instrument will be 
presented in Ravi et al. (in prep.). During the observa- 
tions discussed herein, the DSA-110 was configured as 
described in Ravi et al. (2022). Of particular note is 
that when FRB 20220319D was observed, early in DSA- 
110 commissioning, all candidates at DMs in excess of 
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Figure 1. Dedispersed temporal profile and dynamic spec- 
trum of FRB 20220319D. A value of DMons = 110.95 pc cm? 
was used. The top panel shows the estimated absolute po- 
sition angle based on an approximate parallactic angle, the 
~ 20° error due to RM uncertainty is not included in the 
error bars. The middle panel shows temporal profiles in to- 
tal intensity, linearly polarized intensity, and Stokes V, as la- 
beled. The bottom panel shows the dynamic spectrum of the 
total-intensity data. The time resolution in the time-series 
and dynamic-spectrum plots is 32.768 us, and the temporal 
profile is in units of signal to noise ratio. The reference time 
is MJD 59657.93275696795, which was the burst arrival time 
at OVRO at 1530 MHz. The polarized data on the burst was 
corrected for the measured RM of 50radm7?. 

80 pc cm? were saved for further inspection, regardless 
of the expected DM through the Galactic disk. 

FRB 20220319D was detected during standard com- 
missioning observations, with an arrival time at OVRO 
at 1530 MHz of MJD 59657.93275696795. During these 
observations, the DSA-110 was parked at a pointing- 
center declination of +71.6°. The DSA-110 refer- 
ence position is —118.283? longitude, +37.2334° lat- 
itude. FRB20220319D was detected with a DM of 
110.96 pc cm ?, and a signal to noise ratio (S/N) of 41.7. 
'To derive optimized burst parameters and polarization 
properties, we coherently combined voltage data from 
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the 48 core antennas towards the detection-beam direc- 
tion. The resulting temporal profile and calibrated dy- 
namic spectrum, produced with incoherent dedispersion 
at the native time and frequency resolutions, are shown 
in Figure 1, and optimized parameters are given in Ta- 
ble 1. The polarization analysis was done following pro- 
cedures described in Ravi et al. (2022). The matched- 
filter total-intensity S/N estimate is 79. 

The burst exhibits a narrow temporal profile with 
a strongly modulated spectrum. The DM smearing 
timescale ranges between 8-13 us within the DSA-110 
band, and so the data resolve the temporal structure of 
the burst. Scaling the system-equivalent flux density of 
the DSA-110 by the primary-beam attenuation at the 
burst position, we derive a fluence of 11.2 + 0.2Jy ms. 
Following techniques described in Connor et al. (2020), 
we derive a spectral decorrelation bandwidth, vg = 
3.6 + 0.2 MHz, defined as the 1/e scale of the autocor- 
relation function of the spectrum. The burst spectro- 
temporal morphology is consistent with the population 
of so far non-repeating FRBs (Pleunis et al. 2021). We 
derive a low fractional polarization of just 16% (Fig- 
ure 1), with no significant circular polarization. A Fara- 
day rotation measure (RM) of 50 + 15rad m^? was de- 
tected. The error was estimated by simulating the re- 
covery of RMs with the same linear-polarization S/N as 
FRB 20220319D (Sherman et al., in prep.). 

'The interferometric localization of FRB 20220319D 
was derived following procedures largely described in 
Ravi et al. (2022). Correlation products for all baselines 
formed from the 63 functioning antennas of the DSA-110 
at the time of detection were integrated over 262.1444s 
centered on the burst arrival time. Bandpass calibra- 
tion was done using a 10-min observation of 3C309.1 ob- 
served 11hr prior to the burst detection. Complex gain 
calibration at the time of the burst detection was derived 
using an NRAO VLA Sky Survey (NVSS; Condon et al. 
1998) model and 5 min of visibility data. We show im- 
ages of the point-spread function (PSF; a flat spectrum 
was assumed), and the pre- and post-deconvolution im- 
ages of FRB 20220319D in Figure 2. Briggs weighting 
with a ‘robust’ parameter of 0.5 was used to suppress 
sidelobes, given the high S/N of FRB 20220319D, yield- 
ing a PSF FWHM of 26.9" x 15.6". The position of 
FRB 20220319D, given in Table 1, was derived by fit- 
ting an elliptical Gaussian to the deconvolved image of 
FRB 20220319D. 

We derived the uncertainty in the FRB position us- 
ing data on nine compact bright sources from the Radio 
Fundamental Catalog (RFC; rfc_2022c), obtained within 
2hr of the burst observation. A fit to a linear trend in 
R.A. and Decl. was used to derive arcsecond-level astro- 


Table 1. Properties of FRB 20220319D. 


Parameter Value 
Arrival time (MJD)? 59657.93275696795 
DM (pccm ?) 110.95(1) 

Full-width half-maximum (ms) 0.16(1) 
Fluence (Jy ms) 11.2(2) 

Spectral energy (erg Hz!) 3.3(1) x 1078 
L/1° 0.16(3) 
Iv/1|? 0.04(3) 
RM (rad m?) 50(15) 
va (MHz)^ 3.6(2) 

R.A. (J2000) 02:08:42.7(1) 


Decl. (J2000) +71:02:06.9(6) 


? Arrival time at OVRO at 1530 MHz. 


b EJI is the fraction of linearly polarized fluence, and 
|V/I| is the absolute value of the fraction of circularly 
polarized fluence. 


Cy; is the characteristic spectral decorrelation band- 


width. 


metric corrections, and the corresponding uncertainties. 
These uncertainties (0.52” in R.A. and 0.5” in Decl.) 
were then added in quadrature to the statistical uncer- 
tainty in the fitted burst position (0.29” in R.A. and 
0.17” in Decl.) to derive a final 90% confidence contain- 
ment ellipse with semi-axes of 1.25” in R.A. and 1.18” in 
Decl. This ellipse is shown in the bottom-left panel of 
Figure 2, together with the detrended position offsets of 
the RFC sources. We also assessed the quality of the 
in-field calibration by checking the positions of bright 
(> 50mJy) compact (< 20”) NVSS sources within the 
primary-beam FWHM against NVSS catalog positions. 
These offsets (with the RFC corrections applied) are also 
shown in Figure 2. A more detailed analysis of the DSA- 
110 localization accuracy will be presented in Ravi et al. 
(in prep.). 


3. ASSOCIATION WITH IRAS 02044+7048, AND 
UNCERTAINTIES IN DMism 


The 90% confidence localization region of 
FRB 20220319D is shown in Figure 3 overlaid on a 
Pan-STARRS1 (PS1) i-band image. The galaxy co- 
incident with the FRB location is cataloged in the 
NASA Extragalactic Database as IRAS 02044+7048, 
with a spectroscopic redshift of 0.011 indicating a dis- 
tance of approximately 50Mpc. Using observations 
described below, we derive a spectroscopic redshift of 
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Figure 2. DSA-110 localization of FRB 20220319D. Top left: Point-spread function (PSF) of the DSA-110 for the observation 
of FRB 20220319D, assuming a flat-spectrum source. Top right: Dirty image of FRB 20220319D, with no deconvolution applied. 
Bottom left: Deconvolved image of FRB 20220319D. The synthesized beam is represented by an ellipse with half-power diameters 
of 26.9" x 15.6", at a position angle of 80°. Briggs weighting was used with a ‘robust’ parameter of 0.5, in order to partially 
suppress PSF sidelobes. Baselines shorter than 200A were excluded from the analysis because they were affected by spurious 
cross-coupling. Contours in the preceding three panels are at —0.4, —0.2 (dashed), 0.2, 0.4, 0.6, 0.8 and 0.9 (solid) of the peak 
intensity. The images are centered on the coordinates (R.A. J2000, decl. J2000) = (02:08:42.7, +71:02:06.9). Bottom right: 
Offsets of known astronomical sources from their true positions as measured by the DSA-110. The blue disks show offsets 
of sources from cataloged positions in the NRAO VLA Sky Survey (NVSS; Condon et al. 1998) in a 5 min DSA-110 image 
formed from data taken at a time centered on the burst detection. Only sources with cataloged flux densities > 50mJy and 
measured major-axis diameters < 20", within 2 deg of the pointing center, were considered. The symbol size is proportional to 
flux density. The red crosses show measured offsets of nine RFC calibrator sources (see text for details) observed immediately 
preceding and after the burst. The ellipse indicates the 9096 confidence containment region for FRB 20220319D, derived from 
the RFC calibrators. 
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Figure 3. Pan-STARRS i-band image of IRAS 02044+7048, 
with the 3.7 x 2.2 arcsec 90% confidence localization region 
of FRB 20220319D indicated as a white ellipse. 


0.0111 + 0.0004, indicating a luminosity distance of 
49.6 Mpc. An isophotal fit to the PS1 i-band image in- 
dicates an effective radius of 2.7 + 0.2 kpc, and the FRB 
is offset by 2.3 + 0.3kpc. The isophotal analysis also 
indicates an inclination of 23 + 3deg. A barred spiral 
morphology, with a classification of SBa, is evident from 
the image. 

Under the assumption that FRB 20220319D 
is extragalactic, the spatial association between 
FRB 20220319D and IRAS 02044+7048 is secure. It is 
unlikely that FRB 20220319D lies significantly beyond 
IRAS 02044+7048; the morphology of the galaxy indi- 
cates the presence of a warm-ISM phase, and scattering 
in this ISM would likely cause temporal broadening 
of several hundred milliseconds at our observing fre- 
quencies (Cordes et al. 2022; Ocker et al. 2022). The 
Gravitational Wave Galaxy Catalog (GWGC; White 
et al. 2011) lists 15869 galaxies at distances less than 
50 Mpc, subtending a total of 4deg?. Thus the chance- 
association probability of FRB 20220319D and a galaxy 
at < 50 Mpc is < 1074. 

The extragalactic nature of FRB 20220319D is called 
into question by its DM of 110.95 pccm"?. The pre- 
dicted DM for extragalactic objects along its | = 129.2°, 
b = 49.1? sightline is 132.9pccm~? according to the 
NE2001 model for the Galactic ionized ISM distribu- 
tion (Cordes & Lazio 2002), and 187.7 pc cm ^? accord- 
ing to the YMWI16 model (Yao et al. 2017). The RM 
detection for FRB 20220319D of 5015 rad m7? is, how- 
ever, marginally consistent with being in excess of the 


expected Milky Way contribution along the burst sight- 
line. A recent model for the Galactic RM foreground 
along this sightline (Hutschenreuter et al. 2022) indi- 
cates an expected Galactic RM of —14 + 18rad m^? 
towards FRB 20220319D. We cannot draw conclusions 
from the measured spectral decorrelation bandwidth of 
va = 3.60.2 MHz, which is far in excess of the 0.2 MHz 
decorrelation bandwidth predicted by NE2001 due to 
scattering in the Milky Way ISM, because the origins of 
the spectral structure cannot easily be determined. 

We now analyze the performance of the NE2001 and 
YMW16 models with the aim of determining whether 
the model uncertainties can accommodate an extra- 
galactic origin for FRB 20220319D. Both model the dis- 
tribution of the WIM (Draine 2011) in the Galactic 
disk using observations of radio-wave propagation, with 
a dominant thick-disk component encapsulating several 
overdensities and voids motivated by multiwavelength 
data. The YMW16 model is fit to pulsars with indepen- 
dent DM. and distance estimates, whereas the NE2001 
model is additionally fit to pulsar and extragalactic 
radio-source scattering measurements. For low-latitude 
extragalactic sightlines, the radial extent of the thick 
disk is critical in determining DMīsm. In NE2001, a 
cutoff at a radius of 20 kpc is motivated by observations 
of HII regions in other galaxies and scattering of extra- 
galactic sources. Surveys of Galactic HII regions on the 
other hand motivate the cutoff of 15 kpc in YMWI106. 
Neither cutoff is estimated in the NE2001 or YMW16 
fits; rather they are fixed model choices that do not im- 
pact the match between the models and the data un- 
der consideration. The functional forms for the WIM 
distribution in the thick disk also differ, with NE2001 
incorporating an oblate spheroid, and YMW16 incorpo- 
rating a plane-parallel slab that results in significant de- 
viations from NE2001 at low latitudes in the second and 
third quadrants (e.g., Price et al. 2021). High-latitude 
extragalactic sightlines are likely to be better modeled 
because the scale height of the thick disk is a critical fit- 
ted parameter. However, density fluctuations caused by 
turbulence in the WIM, together with clumps and voids, 
complicate our picture even at high latitudes (Ocker et 
al. 2020). Large-scale WIM inhomogeneities at the disk- 
halo interface such as Galactic chimneys (e.g., Koo et al. 
1992; Normandeau et al. 1996) can result in significant 
departures from spatially smooth models for DMjsm. 

We first note that the recent PSRa sample of pul- 
sar parallax measurements (Deller et al. 2019) reveals 
discrepancies between the NE2001 and YMW16 models 
and measured pulsar distances at low latitudes in the 
second quadrant, where FRB 20220319D appears on the 
sky. In general, known errors in NE2001 and YMW16 
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exhibit significant spatial correlation (e.g., Price et al. 
2021), motivating the consideration of sources along 
sightlines similar to FRB 20220319D. Table 2 lists the 
DMs of the three nearest pulsars in the sample to 
FRB 20220319D, together with the predicted DMs given 
the measured distances. The YMW16 model signifi- 
cantly overpredicts all three DMs, and although NE2001 
is more accurate for the modestly distant pulsars likely 
in the Cygnus-Orion arm, it also overpredicts the DM 
of the more distant PSR J04064-6138. 


A more detailed analysis of the performance of the 
models for sightlines towards distant objects is enabled 
by the recent convergence in the distances to globu- 
lar clusters hosting pulsars. We obtain accurate dis- 
tance estimates (less than few-percent errors) to all 
21 globular clusters with associated pulsars (Manch- 
ester et al. 2005) at |b] > 5? from the compilation 
of Baumgardt & Vasiliev (2021), primarily based on 
Gaia parallaxes (Vasiliev & Baumgardt 2021). Clus- 
ters at lower latitudes are not considered because their 
DMs are influenced by the thin disk and other compo- 
nents in the Galactic plane. For each cluster we com- 
pare the difference between the DMs predicted by the 
NE2001 and YMW16 models for the cluster distances, 
and the true cluster DMs defined as the mean DMs of 
the associated pulsars. The results are shown in Fig- 
ure 4; FRB 20220319D is also included with the pre- 
dicted DM set to DMjsgyy. The errors in NE2001 at 
low latitudes can clearly accommodate an extragalactic 
origin for FRB 20220319D. YMW16 generally performs 
better than NE2001 for the low-latitude clusters, but 
these are nearly all in the first and fourth quadrants 
(with NGC 1851 as the only exception), and so perhaps 
not representative of errors along sightlines closer to the 
Galactic anticenter. 

As an additional check of the Galactic sightline to- 
wards FRB 20220319D, we also consider independent 
measures of DMīsm. First, we use the relation between 
the neutral hydrogen column density and pulsar DMs 
derived by He et al. (2013) to estimate the DMs to- 
wards FRB 20220319D and each globular cluster, as- 
suming no HI gas beyond the clusters. We obtain HI 
column densities from the HI4PI data (HI4PI Collabo- 
ration et al. 2016). The results are shown in Figure 4, 
including errors representing intrinsic scatter in the rela- 
tion. The HI column densities result in underestimates 
of the DMs towards low-latitude clusters, possibly re- 
flecting the increased ionization of HI gas towards the 
central regions of the Galaxy (e.g., Madsen & Reynolds 
2005). The HI-predicted DMjgy for FRB 20220319D 


is consistent with an extragalactic origin, although the 
uncertainties are large. Second, we estimate the DMs 
towards FRB 20220319D and each globular cluster us- 
ing Ha emission measures (EMs) derived from the Wis- 
consin Ha Mapper (WHAM) all sky maps (Haffner et 
al. 2003). We follow the analysis of Berkhuijsen et 
al. (2006), wherein the photon fluxes, Fy, (in units of 
Rayleighs), integrated over the range of Galactic veloci- 
ties are converted to DMs using the following relations: 


EM = 2.25 Fg4 2? 2E(8-V) c7 6 pe (3) 
DM - (EM x D x F). (4) 

f 
F=—— 5 
(1 + €?) (5) 


The expression for F (e.g., Ocker et al. 2020, and 
references therein) encodes an ionized cloudlet model, 
wherein the WIM is composed of cloudlets with a volume 
filling factor f, ¢ = (n?)/(n-)/? ~ 2 captures the cloud 
to cloud variation, and €? € 1 is the density variance 
internal to the cloudlets. Further, E(B — V) is a mea- 
sure of the interstellar redenning along a given sightline, 
and D is the distance through the WIM. We have as- 
sumed a uniform WIM electron temperature of 8000 K, 
and equivalence between the line-of-sight and volume 
filling factors. We find that for the cluster sightlines the 
Ha-based DM estimates perform remarkably well for a 
nominal value F = 0.1 (i.e, f Z 0.4), comparably to 
YMW416 and better than NE2001 for low-latitude sight- 
lines. This is surprising given previous results for pulsar 
DMs (e.g., Schnitzeler 2012), but may be explained by 
globular clusters generally lying beyond the outer extent 
of the WIM, obviating the need for distance corrections 
to the WHAM EMs. For the FRB 20220319D sight- 
line, the Ho flux implies DMīsm = 89 pc cm~’, which is 
22 pccm^? lower than the measured DM. 

We conclude that models for the DM towards 
FRB 20220319D contributed by the Galactic WIM are 
consistent with an extragalactic origin. This se- 
cures the association of the burst with the galaxy 
IRAS 02044+7048. Consistent with the original presen- 
tations of the NE2001 and YMW16 models for the WIM 
distribution, we confirm that there are significant model 
uncertainties at low Galactic latitudes. We used a large 
sample of globular-cluster parallax distances, which have 
DM estimates from associated pulsars, to demonstrate 
that an alternative predictor of DMygm based on the to- 
tal Galactic Ha flux can provide a useful check on the 
aforementioned models. In order to estimate DMīsM 
towards FRB 20220319D we consider pulsars that are 
nearby on the sky and at a similar Galactic latitude. The 
two nearest pulsars within 1 deg of the burst in Galactic 
latitude, PSR J0231+7026 (total separation of 1.9 deg) 
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Table 2. Measured and predicted DMs of PSRz pulsars near FRB 20220319D. 


Pulsar l(deg.) b (deg. FRB offset (deg) Distance (kpc) DM^ NE2001° vMWi16^ 
PSRJ0147--5922 130.1 —2.7 12 2027 i 40.1 31.8 79.0 
PSR J0157+6212 130.6 0.3 9 1.80*0:05 30.2 31.5 59.3 
PSRJ0406--6138 144.0 7.0 15 ASS LES 65.3 123.7 141.7 
? All DMs are given in units of pc cm ?. 
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Figure 4. Difference between the predicted DM and true DM of all globular clusters that host pulsars, sorted by increasing 
|b|. Predictions are made for each globular cluster given its distance and location on the sky, using four different models for the 
warm ISM of the Milky Way. The two literature models are NE2001 (blue solid thick curve) and YMW16 (orange solid thin 
curve). We also convert the measured Ha emission measures (green dashed curve) and HI columns (red dots and shaded region) 
along each sightline to DM predictions using methods described in the text. We also include predictions for FRB 20220319D 


assuming a location beyond the Galactic disk. 


and PSR J0325+6744 (total separation of 7.4 deg), have 
DMs of 46.7 pc cm-? and 65.3 pccm~?, and we consider 
these as conservative and less conservative estimates re- 
spectively of DMysm along the burst sightline. We note 
that the predicted NE2001 and YMW16 DMjsm values 
along the sightlines towards these pulsars are lower than 
along the burst sightline. 


4. THE HOST ENVIRONMENT OF FRB 20220319D 


We obtained long-slit optical spectroscopy of 
IRAS 02044+7048 using the Double Spectrograph 
(DBSP; Oke & Gunn 1982) mounted on the Palo- 
mar 200-inch Hale telescope on 2022 June 01 (UT). 


'The observations were undertaken under conditions of 
1.2" seeing, at an airmass of 1.9, and used a 1.5" slit 
positioned to include both the FRB location and the 
galaxy nucleus at a position angle of 104°. Two 600s 
exposures were obtained with the D55 dichroic and the 
316/7500 grating on the red arm; we only considered 
data from the red arm of the spectrograph given the 
significant extinction along this sightline (Ay = 2.209; 
Schlafly & Finkbeiner 2011). The two-dimensional spec- 
tral data were reduced according to standard procedures 
using a PypeIt-based pipeline (Prochaska et al. 2020). 
We then defined the trace function using observations 
of the standard star Feige 34, and extracted spectra in 
1.5” windows along the slit. Feige34 was also used 
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Table 3. Observed and derived parame- 
ters of the host galaxy of FRB 20220319D, 
IRAS 02044+7048. 1o errors in the last 
significant figures are given in parentheses. 


Parameter Value 

Redshift 0.0111(4) 
Luminosity distance (Mpc) 49.6 
Effective radius (kpc) 2.7(2) 
FRB projected offset (kpc) 2.3(3) 
log M. (Mo) 9.93(7) 

log Z? 0.1(0.2) 

Internal Ay? (utu 

SFR (Mo yr-1)^ 1.8(0.7) 
Inclination (deg) 23(3) 


“Metallicity with respect to solar. 


b V-band extinction corresponding to a 
uniform dust slab. 


© Averaged over the past 100 Myr. 
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Figure 5. Total Ha flux observed at different lo- 


cations within our longslit spectroscopic observations of 
IRAS 02044+7048. Measurements were obtained at 1.5” in- 
tervals along a 1.5” longslit oriented at a position angle of 
104°, which covered both the FRB position and the cen- 
ter of the galaxy. The locations of the galaxy nucleus and 
FRB 20220319D are indicated in the figure. 


for flux calibration, and all spectra were corrected for 
extinction according to the Fitzpatrick & Massa (2007) 
extinction curve. 

We detected Ha emission at several locations along 
the slit, as shown in Figure 5. The strongest Ha emission 
is observed at the galaxy core, from which we also detect 


lines corresponding to [NII] and [SII]. These lines were 
together used to derive a redshift of 0.0111 + 0.0004 for 
IRAS 02044+7048. Although we do not detect H8 or 
[ONI], the ratio log [NII] /Ha = —0.23 + 0.02 indicates a 
< 10% chance of the nuclear ionization being purely due 
to star formation activity (Ho et al. 1997; Kewley et al. 
2006). There is evidence of Ha emission from spiral-arm 
features on either side of the nucleus; these features are 
too extended to represent HII regions. The FRB appears 
associated with the eastern arm, although not at the 
center of its Ha radial profile. No additional evidence 
for or against an arm association can be gleaned from 
the residuals of a fit of the IRAS 02044+7048 image in 
Figure 3 to a two-dimensional Sersic profile. 

In order to determine global properties of 
IRAS 02044+7048, we derived and modeled its spec- 
tral energy distribution (SED). We executed aperture 
photometry on archival images from the Panoramic 
Survey Telescope and Rapid Response System (Pan- 
STARRS; Chambers et al. 2016), Two Micron All Sky 
Survey (2MASS; Skrutskie et al. 2006), and ALLWISE 
(Cutri et al. 2021) surveys. We identified an elliptical 
aperture that captured the i-band extent of the galaxy 
in Pan-STARRS data, and masked likely foreground 
objects. We then used this aperture, convolved with the 
PSF of each survey, to measure extinction-corrected AB 
magnitudes in all filters used by Pan-STARRS, 2MASS, 
and WISE. We modeled this SED using the Prospector 
stellar population synthesis modeling code (Johnson et 
al. 2021). We ran Prospector using recommended tech- 
niques and priors and a non-parametric star-formation 
history (SFH) model (Leja et al. 2019), and sampled 
from the posterior using emcee (Foreman-Mackey et al. 
2013). Non-parametric models for the SFH result in 
less bias in both stellar-mass (M,) and star-formation 
rate (SFR) estimates, because specific SFH models are 
not excluded a priori. We included a model for dust re- 
radiation in the likelihood function. The spectral energy 
distribution of IRAS 020444-7048 is shown in Figure 6, 
together with the results from the Prospector analysis, 
and derived maximum aposteriori probability parame- 
ters are listed in Table 3. The SFR of 1.8 Mc yr^! was 
derived by integrating the SFH over the past 100 Myr. 

IRAS 020444-7048, and the location of 
FRB 20220319D within it, are largely consistent with 
previous results on FRB host galaxies and environ- 
ments (Bhandari et al. 2020; Heintz et al. 2020; Boch- 
enek et al. 2021; Mannings et al. 2021; Bhandari et al. 
2022). Most FRB host galaxies have detectable ongoing 
star formation, and some exhibit signatures of addi- 
tional nuclear ionization sources. The stellar mass of 
IRAS 02044+7048 is entirely consistent with the distri- 
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Figure 6. Left: Prospector fit to the spectral energy distribution (SED) of the host galaxy of FRB 20220319D, 
IRAS 02044--7048. The SED measurements from archival PanSTARRS, 2MASS, and WISE data are shown in red, together 
with representative 10% errors. Filter transmission curves are shown in grey. The maximum aposteriori probability (MAP) 
model photometry is shown in green, together with the MAP spectrum. 100 spectra generated using random draws from the 
posterior distributions of the model parameters are shown in black. Right: Measured SFH of IRAS 02044+7048 (green). The 
results from 100 random draws from the posterior distributions are shown in black. 


bution of masses found for both repeating and so far 
non-repeating FRB hosts. The potential association of 
FRB 20220319D with a spiral-arm feature is consistent 
with the results of Mannings et al. (2021) based on 
Hubble Space Telescope imaging of eight FRB hosts. 
In the absence of very long baseline interferometry, our 
data do not have sufficient angular resolution to deter- 
mine whether or not the FRB originates from within an 
association of young stars (e.g., Tendulkar et al. 2021). 
The location of FRB 20220319D at an offset of ~ 85% of 
the effective radius from the galaxy nucleus is consistent 
with the range of previously observed FRB offsets. 

However, the SFR of IRAS 02044+7048 is rather high 
for its stellar mass, with respect to the observed sam- 
ple of hosts of so far non-repeating FRBs (Bhandari 
et al. 2022). Only one so far non-repeating FRB 
(20190102C, at z = 0.2912) has a comparably high 
specific SFR (sSFR), of logsSFR = —9.74. Although 
IRAS 02044+7048 is consistent with a location on the 
star-forming main sequence of galaxies (Speagle et al. 
2014), non-repeating FRB hosts as a population are con- 
sistent with originating from below this sequence. It 
is possible that FRBs from more actively star-forming 
hosts are selected against in existing radio surveys due to 
propagation effects in the host ISM (Seebeck et al. 2021), 
in which case it is not surprising that IRAS 02044+7048 
has a face-on orientation. 


5. THE MASS OF THE MILKY WAY CGM 


The constraints on the extragalactic DM contribution 
along the sightline towards FRB 20220319D are severe. 
Adopting the two estimates discussed above for DMīsm 
from pulsars nearby to the FRB on the sky, we find 


possible upper limits on the extragalactic DM contribu- 
tion (DMcam +DMigm+DMhost) of either 45.7 pc cm? 
or 64.3pccm-?. Even with our conservative estimates 
for DMjgm, FRB 20220319D has the lowest extragalactic 
DM yet measured for an FRB localized to a host galaxy. 
We thus have the opportunity to stringently bound the 
heretofore poorly constrained value of DMcem, and thus 
directly bound the mass of the Galactic CGM. The IGM 
likely contributes 7 pc cm? towards IRAS 02044+7048, 
assuming an IGM baryon fraction of 0.7 (e.g., Macquart 
et al. 2020). The detection of extended Ha emission 
from the galaxy itself, and in particular at the location 
of FRB 20220319D, indicates that the WIM component 
of the ISM is present, although we do not know where 
within the ISM column the FRB source is located.? We 
simply assume a nominal value of DMpost = 10 pc cm~? 
for our analysis below. Thus, we find that for the 
FRB 20220319D sightline, two possible upper limits on 
DMoaw are 28.7 pccm^? and 47.3 pccm ^ ?, depending 
on the assumption for DMrzsy. These upper limits are 
comparable to previous results from the closest known 
FRB source (FRB 20200120E), in a globular cluster as- 
sociated with M81, which provide possible upper limits 
on DMcem of 32 pc cm? and 42 pc em ?] depending on 
the choice of NE2001 or YMW16 for DMygm (Kirsten et 
al. 2022). The constraints from FRB 20200120E are af- 
fected by a reliance on the NE2001/YMW16 models for 


? With a measurement of scattering in the host ISM, it is possible 
to estimate DMhost (Cordes et al. 2022). However, upper limits 
on the temporal broadening of FRB 20220319D of ~ 0.1 ms are 
not usefully constraining. 
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DMjsm, and by stronger assumptions for the halo DM 
of M81. 

A recent synthesis of models for the Galactic CGM 
DM by Keating & Pen (2020) highlighted the wide range 
of extant predictions. In this work, we use a representa- 
tive set of models to demonstrate the implications of the 
constraint on DMcem from FRB 20220319D. The mod- 
els are illustrated in Figure 7. We also consider measure- 
ments of the DM contributed by sightlines through the 
halo from pulsars in the LMC (Ridley et al. 2013) and 
in the distant (18.5kpc; Baumgardt & Vasiliev 2021) 
globular cluster M53 (Kulkarni et al. 1991; Pan et al. 
2021). Although the LMC is at b ~ —33? and M53 is at 
b ~ +80°, we conservatively subtract the lowest DMrsw 
to be found in either the NE2001 or YMW16 models of 
18 pc cm? to estimate the halo DM contributions along 
their sightlines. In the top two panels of Figure 7, we 
show the resulting upper limits on DMcaw towards M53 
and the LMC, with ranges corresponding to the range of 
associated pulsar DMs. The bottom-left panel shows the 
upper limit on DMcay from FRB 20220319D, placed at 
the virial radius of the Milky Way halo. In all cases we 
show model predictions for DMcgm assuming a spher- 
ically symmetric galactocentric baryon halo, but with 
the DM evaluated along sightlines from the position of 
the Earth (GRAVITY Collaboration et al. 2019). We 
evaluate the modified NFW (mNFW) and Maller & Bul- 
lock (2004) profiles following Prochaska & Zheng (2019), 
with a total Milky Way mass of Miot = 1.5 x 10!2 Mo, 
and a halo baryon fraction of foam = 0.75. The Pen 
(1999) model is evaluated for the same Mis. and with 
the core radius set to the halo virial radius. The semi- 
empirical Miller & Bregman (2013) and Faerman et 
al. (2017) models and the semi-analytic Faerman et al. 
(2020) model are all evaluated as given. Finally, in the 
bottom-right panel of Figure 7, we also show the model 
predictions for the halo baryon density along the LMC 
sightline, together with constraints from ram-pressure 
stripping analyses of dwarf galaxies (Salem et al. 2015; 
Putman et al. 2021) and modeling of OVII and OVIII 
emission measurements (Kaaret et al. 2020). 

'The data are in favor of models that predict lower to- 
tal CGM DMs and lower inner densities (e.g., Pen 1999; 
Faerman et al. 2020; Miller & Bregman 2013). The DMs 
of M53 pulsars and FRB 20220319D deliver consistent 
constraints on the range of models, albeit on very dif- 
ferent radial distance scales. The constraints from LMC 
pulsars are less impactful, but nonetheless also exclude 
the three models (Maller & Bullock 2004; Faerman et 
al. 2017; Prochaska & Zheng 2019) considered here that 
predict larger values of DMcem. We note that our treat- 
ment of the CGM contribution to the DMs of LMC pul- 


sars is far more conservative than that of Anderson & 
Bregman (2010); this is motivated by the uncertainties 
discussed above in estimating DMrsy. The DM con- 
straints are consistent with the most robust density es- 
timates that we consider (Salem et al. 2015; Kaaret et 
al. 2020). The density estimates are however subject 
to the uncertainties discussed in 81, and we proceed by 
considering only the constraints on DMcay. 

We now derive constraints on the mass of the Galac- 
tic CGM, Mcem, by finding 8 models that are consis- 
tent with the FRB 20220319D and M53 DMs. When 
converting the electron column density to a mass col- 
umn density, we assume 1.18 proton masses per elec- 
tron following Yamasaki & Totani (2020), which corre- 
sponds to roughly solar metallicity gas. The £ profile is 
widely adopted in the field to convert (column-)density 
estimates to halo masses (e.g., Miller & Bregman 2013; 
Salem et al. 2015; Kaaret et al. 2020), as it is empirically 
motivated. The electron number density at a radius r is 
given by 


EV —3B/2 

uc | ED. 
where ng is a central density, and r, is a core radius that 
we fix to 0.47 kpc following Miller & Bregman (2013). 
The exact value of the core radius is unimportant for 
our conclusions regarding Mcem. ‘Typical values of 
B found using various CGM tracers are in the range 
0.4 S B S 0.5 (e.g., Miller & Bregman 2013, 2015; Salem 
et al. 2015; Kaaret et al. 2020). We consider three con- 
straints on DMcaw: the lower and upper constraints 
from FRB 20220319D, and the lower constraint from the 
M53 pulsars. For each constraint, we derive correspond- 
ing values of Mccgm from an integrated 8 profile for dif- 
ferent values of 8. A useful joint constraint on Meem 
and f is not possible with the data in hand. The results 
are shown in Figure 8. 

For typical values of Ø of between 0.4 and 0.5, 
the DM of FRB 20220319D implies an upper limit on 
log Mcam of between 11.0 and 10.8. The M53 pul- 
sars are even more constraining, limiting log Meam to 
below a value of between 10.7 to 10.5. We consider 
the FRB 20220319D constraints to be more conserva- 
tive because the M53 estimate of DMcam is very sen- 
sitive to the assumed DMjgm along its sightline. The 
mass constraints all favor values of foam < 0.5 (re- 
call foam = Meeutiac) for our fiducial value of Mito: = 


n(r) = no 


1.5x 10? Mo, with M53 potentially accommodating val- 
ues of foam S 0.2 for reasonable values of 5. 
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Figure 7. Illustration of the constraining power of different measurements of the CGM DM and density on a range of models. 
In all panels, curves show predictions of the default modified NFW (mNFW) model of Prochaska & Zheng (2019) (blue), the 
model of Maller & Bullock (2004) (orange; MB04), the model of Miller & Bregman (2013) (green; MB13), a Pen (1999) model 
(Pen99) assuming a core radius set to the virial radius (red), and the models of Faerman et al. (2017) (purple; Faermanl) and 
Faerman et al. (2020) (brown; Faerman2). Top panels: upper limit on the CGM DM towards M53 (left) and the LMC (right). 
In both cases, the minimum DM (18 pc cem ?) across the sky in either the NE2001 or YMW16 models has been subtracted. The 
ranges indicate the range of DMs of pulsars associated with each object. Bottom left: upper limit on the CGM DM towards 
FRB 20220319D. The range indicates different assumptions for the ISM DM that correspond to the two closest pulsars on the 
sky. Bottom right: constraints on the CGM density compared with the models. The dark disks indicate measurements using 
ram-pressure stripping modeling of the LMC (Salem et al. 2015) and modeling of diffuse X-ray line emission (Kaaret et al. 
2020). The lower limits displayed as triangles indicate estimates using ram-pressure stripping modeling of a large selection of 


Milky Way satellites (Putman et al. 2021). 
6. DISCUSSION 


Our results demonstrate that the total Galactic 
baryon mass is likely significantly lower than the cos- 
mological average for a halo as massive as the Milky 
Way. The mass constraints we derive are consistent 
with some previous observational results from analyses 
of OVII and OVIII emission and absorption (e.g., Li 
& Bregman 2017; Kaaret et al. 2020) and ram-pressure 
stripping of the LMC (Salem et al. 2015), but inconsis- 
tent with other results that posit higher CGM masses 
(e.g., Faerman et al. 2017; Yamasaki & Totani 2020). 
Analyses of the dynamics of Milky Way satellites im- 
ply a total virial mass of Mig, = (1.4 0.3) x 101? Mo 
(Watkins et al. 2010), a recent Gaia-based analysis 
of the dynamics of Milky Way globular clusters yields 
My = (1.32: 0.3) x 101? Mo (Posti & Helmi 2019), and 
an analysis of the dynamics of the Magellanic Stream 
yields Mo, = (1.5 + 0.3) x 101? M (Craig et al. 2022). 
It is possible that a full accounting for the orbit of the 


LMC would reduce the estimates of Mio, by ~ 15% 
(Correa Magnus & Vasiliev 2022). However, in all 
cases foam S 0.4 is implied by FRB 20220319D, and 
foam S 0.2 is implied by the M53 DM. Our constraints 
are affected by different systematic effects to previous 
results, and do not suffer from uncertainties in model- 
ing the chemical and thermal states of the CGM, nor 
from uncertainties in modeling the interaction of satel- 
lite galaxies with the CGM. These values are consistent 
with several simulations of the CGM contents of galaxies 
like the Milky Way that account for the effects of kinetic 
and thermal feedback on reducing the halo baryon con- 
tent. 

The observational constraints on Mcam are affected 
by uncertainties in identifying the CGM DM contri- 
bution along the sightlines of interest. For the low- 
latitude FRB 20220319D sightline, we adopted conser- 
vative estimates of DMygm, and there are uncertainties 
at the level of a few pccm-? in the DMjaw and DMyost 
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Figure 8. Upper limits on the CGM mass for different 
assumptions for the radial density profile. We adopt a 8 
model (see text for details), and show limits corresponding 
to the upper and lower ends of the range of constraints from 
FRB 20220319D (dashed and solid lines respectively), and 
from the most constraining pulsar in M53 (red solid line). 
Typical values of 6 found in the literature are = 0.4. For 
these values, the inferred foam values are typically below 
0.5, as indicated by blue horizontal lines. 


terms. We therefore consider the constraints on Mcam 
based on FRB 20220319D to be robust to uncertainties 
in DMcem. The impact of uncertainty in DMrsm for the 
M53 sightline is greater given the low values of DMcam 
under consideration. Although we attempt to be con- 
servative, even a few-unit increase in DMcgm would sig- 
nificantly raise the derived upper limits on Meem. Our 
analysis is also sensitive to how we model the distribu- 
tion of baryons in the CGM. First, we do not account 
for sightline to sightline variance, which simulations sug- 
gest may be at the level of = 10% (e.g., Zheng et al. 
2020; Ramesh et al. 2022). Second, it may be that a 6 
profile is not the correct model to adopt, and our con- 
straints are sensitive to the specific value of 3. Guidance 
on these points from simulations for future analyses of 
DMcem will be important. Finally, although we assume 
a specific helium abundance in translating the DM con- 
straints to a total CGM mass, our results are sensitive 
to variations in this assumption only at the < 10% level. 

The discovery of FRB 20220319D has major implica- 
tions for the search for FRBs in the local universe, and 
for future FRB studies of the CGM of the Milky Way 
and other galaxies. First, we have shown that for low 
Galactic latitudes (|b| S; 15°; see Figure 4) there can 
be significant uncertainties in the widely used NE2001 
and YMW16 models for DMīsm, which may lead to 
nearby extragalactic FRBs being missed in surveys. In 
all cases we recommend that FRB surveys save candi- 
dates at DMs below the model values of DMyjsy. The 
possibility of missing FRBs with low extragalactic DMs 


also needs to be considered when using FRB samples to 
directly infer a characteristic DMcem (e.g., Platts et al. 
2020). These studies are also likely to be impacted by 
the uncertainties in deriving DMcem that in most cases 
will be on the order of the values of DMcam allowed 
by our analysis (see Figure 7). We recommend that in- 
ferences of DMcam with FRBs focus on high-latitude 
sightlines towards FRBs interferometrically localized to 
nearby galaxies, wherein uncertainties in DMjgm and 
DMhnost can be minimized. Nonetheless, if the charac- 
teristic DMcam is indeed on the order of 10pccm-?, 
a very large number of local FRBs will be required to 
suppress variance between halo sightlines, and uncer- 
tainties in other DM contributions. These samples may 
be furnished by future coherent all-sky radio monitors 
(e.g., Lin et al. 2022). Finally, consistent with previ- 
ous studies (Keating & Pen 2020; Kirsten et al. 2022), 
our constraints on foam firmly exclude the fiducial pa- 
rameterizations of the widely used modified-NFW and 
Maller & Bullock (2004) models for the CGM DM as 
described by Prochaska & Zheng (2019). This impacts 
studies that have assumed the correspondingly large val- 
ues of DMcem (e.g., Ravi 2019), as well as forecasts for 
the contributions of the CGM of intervening galaxies to 
FRB DMs. 

The host of FRB 20220319D, IRAS 02044+7048, ap- 
pears to have an unusually large specific SFR. This is 
despite FRB 20220319D being the closest so far non- 
repeating FRB yet discovered, and the strong evolution 
towards higher SFRs of the star-forming main sequence 
of galaxies with redshift (Speagle et al. 2014). Our anal- 
ysis of the SFH of the IRAS 02044+7048 host galaxy 
suggests that most stars were formed in a burst in the 
last S 2Gyr. Thus it is likely that the progenitor of 
FRB 20220319D does not have a delay time (i.e., age 
from its formation epoch) in excess of this timescale. 
Although several progenitor scenarios remain consistent 
with this constraint, the additional possible association 
of the FRB source with a spiral arm may imply a pro- 
genitor age that is a fraction of a ~ 100 Myr dynamical 
timescale. We note that it is not straightforward to com- 
pare our inferred SFR, which is averaged over the past 
100 Myr of the SFH, to the Ha-based SFRs common in 
the FRB literature (e.g., Bhandari et al. 2022), which 
are sensitive to the last ~ 10 Myr (Kennicutt 1998). A 
direct comparison of non-parametric SFH estimates may 
lead to a more complete view of the formation channels 
of FRB progenitors. 


7. SUMMARY AND CONCLUSIONS 


We present the DSA-110 discovery and interfero- 
metric localization of the nearest so far non-repeating 
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FRB (20220319). The burst was observed with a high 
S/N of 79, a DM of 110.95 pccm ?, and has a lin- 
ear polarization fraction of 16 + 3% and an RM of 
50 + 15radm ? (Figure 1, and Table 1). The FRB 
originated from the position (R.A. J2000, decl. J2000) 
= (02:08:42.7(1), +71:02:06.9(6)), where uncertainties 
in the last significant figures are given in parenthe- 
ses (Figure 2). We associate FRB 20220319D with a 
spiral-arm feature of a face-on star-forming galaxy at a 
distance of 50 Mpc (IRAS 02044+7048; Figures 3 and 
5; Table 3). IRAS020444-7048 is moderately mas- 
sive (log M. = 9.93 + 0.07), with an approximately 
solar metallicity and an SFR. averaged over the past 
100 Myr of 1.8 + 0.7 Mo yr ! (Figure 6). The low DM 
of FRB 20220319D motivates a sensitive new constraint 
on the baryon mass of the Milky Way CGM. Our con- 
clusions are summarized as follows. 


e The spatial coincidence of FRB 20220319D and 
IRAS020444-7048 (false alarm probability of 
< 10-^), and a potential RM detection for 
FRB20220319D that is in excess of the Galac- 
tic foreground, both imply a secure association. 
However, the DM of FRB 20220319D is somewhat 
lower than the predicted Galactic ISM DM along 
its sightline from leading models for the WIM dis- 
tribution (NE2001 and YMW16; Cordes & Lazio 
2002; Yao et al. 2017). We show through an anal- 
ysis of the DMs towards pulsar-hosting globular 
clusters with Gaia parallax distances (Figure 4) 
that there are significant uncertainties in NE2001 
predictions for DMysm for |b| < 15°, and that 
YMW41G6 is particularly uncertain at low latitudes 
in the second quadrant, where FRB 20220319D is 
located (see also Table 2). We find that estimates 
for DMygsm based on Ha EMs perform compara- 
bly well to the models. We conclude that models 
for DMīsm towards FRB 20220319D can accom- 
modate an extragalactic origin, and urge FRB sur- 
veys to consider candidate bursts at DMs below 
standard model predictions for DMjsm. 


The SFH of IRAS 02044+7048 is consistent with 
most stars being formed in a burst within the 
past two Gigayears. This, together with the pos- 
sible location of FRB 20220319D within a spi- 
ral arm of IRAS 02044+7048, suggests a mod- 
erately low delay time for the burst progenitor. 
Despite its nearby distance, the specific SFR of 
IRAS 02044+7048 derived from SED modeling is 
larger than all but one so far non-repeating FRB, 
although the comparison with literature samples 


of host galaxies would be aided by the widespread 
use of non-parametric SFH analyses. 


We find conservative upper limits on the Galactic 
CGM contribution to the DM of FRB 20220319D 
of either 28.7 pc cm? or 47.3 pc em ?: the two val- 
ues are based on DMysy-estimates from the DMs 
of two pulsars nearby on the sky. These limits ex- 
clude some literature models (Maller & Bullock 
2004; Faerman et al. 2017; Prochaska & Zheng 
2019) for the baryon distribution in the CGM (Fig- 
ure 7). We find that consistent results are ob- 
tained from estimates of the CGM DM contribu- 
tions towards pulsars in the LMC and in the dis- 
tant (18.5 kpc) high-latitude globular cluster M53. 
The M53 pulsars in particular provide interest- 
ing constraints in the inner regions of the Milky 
Way halo that complement the FRB constraint, 
although the M53 DMcem estimate is sensitive to 
the assumed ISM contribution. 


We derive an upper limit on the mass of baryons in 
the Galactic CGM of log Mcem S 10.8 — 11 from 
FRB 20220319D, and log Mcem S 10.5— 10.7 from 
M53 (Figure 8). The range corresponds to the 
range of assumed indices (0.4-0.5) of the 8 pro- 
file for the baryon halo density distribution. For 
a fiducial total mass (baryons and dark matter) 
of the Milky Way of 1.5 x 10? M5, and an as- 
sumed baryon-disk mass of 6 x 101? M, the Milky 
Way contains « 6096 of the cosmological-average 
baryon mass. This is consistent with a baryon 
census, predicted by galaxy-formation simulations, 
wherein these “missing” baryons are expelled from 
halos like the Milky Way into the IGM through ki- 
netic and thermal feedback from AGN, supernovae 
and massive stars. Although our analysis relies 
on just a few sightlines, the conservative upper 
bounds on Mcam assuming spherical symmetry 
are robust to variations in the density distribution 
of halo baryons, and to assumptions on the chem- 
ical and thermal state of the halo. 


Studies of unseen baryons with FRB propagation sig- 
natures, like DMs, promise to transform our understand- 
ing of the distribution of baryons around and in be- 
tween galaxies (Ravi et al. 2019). As more FRBs like 
FRB 20220319D are localized to nearby galaxies (see 
also Kirsten et al. 2022), constraints on the CGM con- 
tent of the Milky Way will continue to improve, and 
direct measurements of Mcaw may be possible. These 
measurements are likely to be important in assembling 
an in-situ picture of the processes whereby a galaxy 
grows from and impacts its environment. 
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